Second-harmonic generation and linear electro-optical coefficients of SiC polytypes 

and nanotubes 
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The second-order nonlinear optical susceptibility (xitc) anc ^ unear electro-optical coefficient (r a & c ) 
of a large number of single-walled zigzag, armchair and chiral SiC nanotubes (SiC-NTs) as well as 
bulk SiC polytypes (2H-, 4H-, 6H- and 3C-SiC) and single graphitic SiC sheet have been calculated 
from first-principles. The calculations are based on density functional theory in the local density 
approximation and highly accurate full-potential projector augmented-wave method is used. Both 
the zigzag and chiral SiC-NTs are found to exhibit large second-order nonlinear optical behavior with 

(2) 

the Xabc an d r abc coefficients being up to ten-times larger than that of bulk SiC polytypes, and also 
being up to thirteen-times larger than the counterparts of the corresponding BN-NTs, indicating 
that SiC-NTs are promising materials for nonlinear optical and opto-electric applications. The 
prominant features in the spectra of ^i? (— 2w, w, u) of the SiC-NTs are correlated with the features 
in the linear optical dielectric function e{ui) in terms of single-photon and two-photon resonances. 

PACS numbers: 42.65.-k,42.70.Nq,73.22.-f,78.67.Ch 



I. INTRODUCTION 

Carbon nanotubes (CNTs) have attracted considerable 
interest worldwide ever since their discovery in 1991 [l|, 
mainly because of their unusual properties and great po- 
tentials for technological applications. CNTs can be re- 
garded as a layer of graphene sheet rolled up in a tubular 
form 0], and the structure of a CNT is specified by a 
chiral vector defined by a pair of integers (n, to). CNTs 
can be chiral or nonchiral depending on the way they are 
rolled up. Their physical properties, in particular, optical 
dielectric functions, depend sensitively on their chirality, 
i.e., the (n,m) indices (see, e.g., Refs. land! and ref- 
erences therein). Apart from CNTs, inorganic tubular 
materials, such as BN A1N Q, GaN [|], have also 

been predicted and synthesized. These tubular materials 
also display some very interesting properties distinctly 
different from their bulks. 

Bulk silicon carbide (SiC) crystallizes in either a cubic 
or a hexagonal form, and exbihits polytypism IToj |. 
The polytypes are made of identical hexagonal layers 
with different stacking sequences. These polytypes are 
semiconductors with a range of band gaps, from 2.39 
cV in the zincblende (3C) to 3.33 eV in the wurtzite 
polytype (2H) [|, [Tl|. Furthermore, 3C- and 6H-SiC 
are used for high temperature, high-power and high- 
frequency devices [H, IH, [T3, [l5| due to their unique 
properties while 6H-SiC with a band gap of 2.86 
eV is a useful material for blue light-emitting diode ap- 
plications [ll[. Recently, SiC-NTs were also synthesized 
by the reaction of CNTs and SiO at different tempera- 
tures |l7j. This has stimulated a number of theoretical 
and experimental investigations on the tubular form of 
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the SiC [ialial2a|2ll|22|. Based on density-functional 
calculations, Miyamoto and Yu *23| predicted the ex- 
istence of graphitic and tubular forms of SiC and also 
proposed their synthesis using an extreme hole injection 
technique. They also reported that the strain energies of 
SiC-NTs are lower than that of CNTs, and that the band 
gaps of SiC-NTs can be direct or indirect, depending 
on the chirality. Using both tight-binding molecular dy- 
namics and ab initio methods, Menon and coworkers [201 ] 
showed that single- walled SiC-NTs are highly stable with 
a large band gap. Zhao et al. [24[ also investigated the- 
oretically the strain energy, atomic and electronic struc- 
ture of SiC-NTs with or without hydrogenation. Gali 
performed an ab initio study of the effect of nitrogen 
and boron impurities on the band structure of the SiC- 
NTs 

Unlike CNTs, SiC-NTs are polar materials and there- 
fore, may exhibit some unusual physical properties that 
CNTs may not have. For example, like BN-NTs, [H ^ 
zigzag SiC-NTs may become piezoelectric, and also show 
second-order non-linear optical response. Despite the in- 
tensive theoretical studies mentioned above, no ab initio 
calculation of the dielectric response and optical prop- 
erties of SiC-NTs has been reported, perhaps because 
of the heavy demand on computer resources. A knowl- 
edge of the optical properties of SiC-NTs is important for 
their optical and electrooptical applications. Therefore, 
we have recently carried out a series of ab initio calcula- 
tions in order to analyze the linear optical features and 
underlying band structure of all three types of the SiC- 
NTs as well as their possible dependence on diameter and 
chirality. [28] In this work, we investigate the second-order 
optical susceptibility and also linear electro-optical coef- 
ficient of the SiC-NTs as well as bulk SiC polytypes. The 
primary objective of this work is to find out the features 
and magnitude of the second-harmonic generation and 
linear electro-optical coefficients of the SiC-NTs in or- 
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der to see whether they have any potential applications 
in nonlinear optical and electro-optical devices such as 
second-harmonic generation, sum-frequency generation, 
and electrical optical switch. The second objective is 
to identify distinguished differences in nonlinear optical 
properties between SiC-NTs and CNTs and also be- 
tween SiC-NTs and BN-NTs [13 . 

The rest of this paper is organized as follows. In Sec. 
II, the theoretical approach and computational method 
are described. In Sec. Ill, the calculated second-order 
nonlinear optical susceptibility and linear electro-optical 
coefficients of SiC nanotubes are presented. For com- 
parison, the second-order nonlinear optical susceptibility 
and linear electro-optical coefficients of bulk SiC poly- 
types and single graphitic SiC sheet are also calculated 
and presented. Finally, in Sec. IV, a summary is given. 



II. THEORY AND COMPUTATIONAL 
METHOD 



Our ab initio calculations for the SiC-NTs were per- 
formed using highly accurate full-potential projector 
augmented- wave (PAW) method [29| . as implemented in 
the vasp package [3(J. They are based on density func- 
tional theory (DFT) with the local density approxima- 
tion (LDA). A supercell geometry was adopted so that 
the nanotubes are aligned in a square array with the clos- 
est distance between adjacent nanotubes being at least 
10 A. A large plane- wave energy cutoff of 450 eV was 
used throughout. We consider a large number of rep- 
resentative SiC-NTs with a range of diameters from all 
three types, as listed in Ref. [H. 

First, the ideal nanotubes were constructed by rolling- 
up a hexagonal SiC sheet. The atomic positions and 
lattice constants were then fully relaxed by the conjugate 
gradient technique. The theoretical equilibrium atomic 
positions and lattice constant were obtained when the 
forces acting on all the atoms and the uniaxial stress 
were less than 0.04 eV/A and 1.0 kBar, respectively. The 
calculated equilibrium lattice constants T and curvature 
energies E c (total energy relative to that of single SiC 
sheet) as well as the other computational details have 
been reported before [28| . 

The final self-consistent electronic band structure cal- 
culations were then carried out for the theoretically de- 
termined SiC-NT structures. In this work, the non- 
linear optical properties were calculated based on the 
independent-particle approximation, i.e., the excitonic 
effects and the local-field corrections were neglected. 
Note that the band gap of SiC in the wurtzite (2H) 
structure from our LDA calculations is 2.12 eV (see Sec. 
IIIA below), being 1.18 eV smaller than the experimen- 
tal value [1, [ll|. This suggests that the quasi-particle 
self-energy corrections to the optical peak positions may 
amount to ~1 cV in the SiC systems. Furthermore, 
because the LDA underestimates the energy gaps, the 
calculated dielectric function and second-order nonlinear 



optical susceptibility in the static limit might be slightly 
too large in general. For example, our calculated dielec- 
tric function at eV for the electric field parallel and 
perpendicular to the c-axis of 2H-SiC is 8.33 and 7.90, 
respectively, being slightly larger than the corresponding 
measured eoo, 6.84 and 6.51. @ Therefore, the results 
from the present LDA-independent-particle calculations 
might not be able to be compared quantitatively with 
experiments, though they would certainly be useful to 
study the trends and characteristics of the optical prop- 
erties of the SiC-NTs. Nonetheless, our previous calcu- 
lations show that the dielectric functions of graphite [|[ 
and also of h-BN [27J calculated within the independent- 
particle picture are in reasonable agreement with experi- 
ments, perhaps because of the accidental cancellation of 
the self-energy corrections by the excitonic effects. 

Following previous nonlinear optical calculations 3, 
1271 [3l| . the imaginary part of the second-order opti- 
cal susceptibility due to direction interband transitons 
is given by [UJ 
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where the contribution due to the so-called virtual- 
electron (VE) process is 
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dipole transition matrix elements pfj —< kj|p a |kz > were 
obtained from the self-consistent band structures within 
the PAW formalism [33| . The real part of the second- 
order optical susceptibility is then obtained from x'alc 
by a Kramer-Kronig transformation 

'^^y^. ( 4) 



X ' (2) (-2w,w,w) 



The linear electro-optic coefficient r a b c (uj) is connected 

(2) 

to the second-order optical susceptibility x„(, c (-w, uj, 0) 
through the relation [34| 

xite(-w>w,0) = --nl(w)nl(w)r abc (w) (5) 

where n a (uj) is the refraction index in the a-direction. In 
the zero frequency limit, 



lim n Xabc(- 2uj i w > w ) = lim n Xo6c(~ W ) W > °)- ( 6 ) 
uj — »0 to — ►U 
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Therefore, 

rabc(0) = " ^(oKto) ^^'" 2 ^'^' (7) 

Furthermore, for the photon energy ftw well below the 
band gap, the linear electro-optic coefficient r a b c (uj) ~ 

( 2) 

^afcc(O) because X a b c (~ ^oj, u>, uj) and n(cu) are nearly con- 
stant in this low fre quen cy region, as shown in the next 
section and in Ref. [281 ] . 

In the present calculations, the ^-function in Eqs. (2)- 
(3) is approximated by a Gaussian function with T = 0.2 
eV. A uniform fc-point grid (Ixlxn) along the nan- 
otube axis (z-axis) with the number n being from 40 

'(2) 

to 100, is used. Furthermore, to ensure that Xabc ca ^ 
culated via Kramer-Kronig transformation (Eq. (4)) is 
reliable, at least ten energy bands per atom are included 
in the present optical calculations. The unit cell vol- 
ume fi in Eqs. (2)-(3) is not well defined for nanotubes. 
Therefore, like the previous calculations [!, d, H3, l28l ]. 
we used the effective unit cell volume of the nanotubes 
rather than the volume of the supercells which is arbi- 
trary. The effective unit cell of a nanotube is given by 
fl = n[(D/2 + d/2) 2 - {D/2 - d/2) 2 ]T = nDdT where 
d is the thickness of the nanotube cylinder which is set 
to the interlayer distance of /i-SiC (3.51 A [Hf). D and 
T are the diameter and length of translational vector of 
the nanotube [28|], respectively. 

III. RESULTS AND DISCUSSION 
A. Bulk SiC polytypes 

For comparison with the SiC-NTs, we first investigated 
the second-order nonlinear optical susceptibility of bulk 
SiC polytypes. SiC exists in either a cubic or a hexagonal 
polymorphic structure. The polytypes are made of iden- 
tical hexagonal layers with different stacking sequences. 
SiC occurs in nearly 200 polytypes. (3f| Here we consider 
only the four common 2H, 4H, 6H and 3C types of SiC. 
The zincblende (3C) structure, with pure cubic stacking 
of the Si-C double layers in the [111] direction, is one 
of the two most extreme polytypes. The other is the 
wurtzite (2H) structure with pure hexagonal stacking in 
the [0001] direction. The atomic positions and lattice 
constants of the four polytypes were theoretically opti- 
mized. In the atomic structure optimizations, a uniform 
fc-point grid of 30 x 30 x n with n being from 15 to 30 was 
used. The calculated equilibrium lattice constants of the 
SiC polytypes are listed in Table I, together with previous 
theoretical and experimental values. Table I shows that 
our theoretical lattice constants are in very good agree- 
ment with the previous calculations [H| , though they are 
slightly smaller (within 0.7 %) than the corresponding 
experimental values [9']. Our calculated indirect band 
gaps for the 3C, 2H, 4H and 6H types of SiC are, re- 
spectively, 1.29 eV, 2.12 eV, 2.17 eV and 1.96 eV. These 



band gap values are smaller by about 1.1 eV than the 
corresponding experimental values @ of 2.39 (3C), 3.33 
eV (2H), 3.26 eV (4H), and 2.86 eV (6H). Nonetheless, 
our theoretical band gap values are in very good agree- 
ment with previous LDA calculations (see, e.g., Ref. [13 
and references therein). The difference in the theoretical 
DFT-LDA and experimental band gap values can be at- 
tributed to the quasi-particle self-energy corrections, as 
has been shown explicitly in previous quasiparticle band 
structure calculations for the SiC polytypes [37| . 

In Fig. 1, we display the imaginary and real parts 
of three nonvanishing components of the second-order 
nonlinear optical susceptibility (tensor) \^ ( — uj, uj) 
of the three hexagonal SiC polytypes. The imagi- 
nary and real parts of the nonvanishing component 

(2) 

Xxyz(—2u>,u>,uj) of the cubic SiC are plotted in Fig. 2. 
In the nonlinear optical calculations, we used a denser fc- 
point grid of 50 x 50 x n with n being from 10 to 30 for the 
three hexagonal SiC polytypes, and n = 50 for the cubic 

(2) 

SiC. Importantly, we find numerically that only the Xzzz, 

(2) (2) (2) (2) (2) , (2) , 

Xzxx, Xzyy, Xxzx, Xyzy, Xxxz, and Xyyz components are 
nonzero for the three hexagonal polytypes. Moreoever, 

(2) (2) i (2) (2) (2) (2) rr U ■ ■ 

Xzyy = Xzxx, and Xyyz = Xxxz = Xxzx = XyAt, lnls 1S 

consistent with the symmetry consideration, [381 ] demon- 
strating that our numerical method and calculations are 

(2) (2) 

qualitatively correct. For the 3C-SiC, only Xxyz, Xzxy, 

(2) 

Xyzx components are nonzero, and all the three com- 
ponents are numerically found to be identical, as they 
should. Fig. 1 shows that the line shape, amplitude and 

(2) (2) 

feature positions of the Xzxx and Xxzx components for 
all the hexagonal SiC polytypes are very similar. On the 
other hand, the Xzzz component varies significantly as 
one moves from 2H- to 6H-SiC. This is mainly because 
in this case all the electric fields are polarized parallel 

(2\ 

to the z-axis and hence the Xzzz component should be 
rather sensitive to the layer stacking sequence along the 

(2) 

z-axis. We also note that the shape of the Xx% U of 3C SiC 

(2) 

is similar to that of the Xzzz of 6H-SiC (Fig. 1 and Fig. 

(2) 

2). Nonetheless, the features in the Xxyz of the cubic 
structure appears to have the largest magnitude among 
the four SiC polytypes studied here. 

We find that our calculated nonzero components of 
•^( 2 )(_ 2ui,ui,ui) of all the four polytypes agree quite well 
with the recent DFT-LDA calculations by Adolph and 
Bechstedt [39|. For detailed comparison, we plot the 
imaginary part of the nonzero components of SiC in the 
cubic and wurtzite structures from the previous calcula- 
tions by Adolph and Bechstedt [1^] in Fig. 2 (cubic) and 
Fig. 3 (wurtzite), together with our results. One can 

(2) 

see from Figs. 2 and 3(c) that Xxyz(— 2w, ui, w) of 3C- 

(2) 

SiC and Xxzx(— 2ui, w, uj) of 2H-SiC from the present and 
previous calculations are almost identical. Nonetheless, 

(2\ 

there arc some discernible discrepancies in the \zxx 

and 

(2) 

Xzzz components of 2H-SiC between the two calculations 
[Fig. 3(a-b)]. 
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TABLE I: Calculated lattice constants of SiC polytypes, together with other theoretical (Cal.) (Bechstedt et al. [36J) and 
experimental (Exp.) (Ref. [Sj) lattice constants. 





2H 


4H 


6H 


3C 




a (A), c (A) 


a (A), c (A) 


a (A), c (A) 


a (A) 


This work 


3.057, 5.018 


3.060, 10.018 


3.061, 15.018 


4.331 


Cal. 


3.057, 5.016 


3.061, 10.012 


3.062, 15.012 


4.332 


Exp. 


3.076, 5.048 


3.073, 10.053 


3.081, 15.117 


4.360 
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FIG. 1: (Color online) Real and imaginary parts of \^ (— 2uj, ui, ui) of the hexagonal SiC polytypes. 
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FIG. 2: (Color online) Real and imaginary parts of 
Xxyz(—2iu,uj, lo) of SiC in the zinc-blende (3C) structure. Pre- 
vious calculations by Adolph et al. [3^] are also displayed for 
comparison. 



In Table II, the calculated zero frequency linear electro- 
optic coefficient r(0) as well as the corresponding second- 



TABLE II: Calculated static refraction index n, second- 
order optical susceptibility x (0) (pm/V) and linear electro- 
optical coefficient r a bc (pm/V) of the hexagonal SiC polytypes 
(a), zinc- blende SiC (b) and also isolated SiC sheet (c). 



(a) 


n x (n z ) 


v (2) y (2) 


fzzz-) Tzxx 


2H-S1C 


2.81 (2.88) 


11.5, -17.9 


-0.33, 0.55 


4H-SiC 


2.83 (2.87) 


32.8, -20.4 


-0.97, 0.62 


6H-SiC 


2.82 (2.89) 


38.6, -20.4 


-1.11, 0.61 


(b) 


n 


(2) 
Xxyz 


Txyz 


3C-S1C 


2.95 


34.2 


-0.90 


(c) 


n a (n c ) 


(2) (2) 
v v 

A.bbb~> A-aab 


Tbbbi Taab 


SiC sheet 


3.19 (2.01) 


227.8, -227.8 


-11.08, 11.08 



order nonlinear optical susceptibility x (0, 0) 0) are 
listed. The r(0) is calculated from the corresponding 
X (2) (0,0,0) by using Eq. (7). We find that the calcu- 
lated xizx(0) and xiSr(0) f° r the three hexagonal SiC 
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polytypes differs slightly from each other, i.e., only ap- 
proximately satisfying the requirement of the Kleiman 

symmetry In Table II, therefore, the xixx(O) is listed 

(2) (2) 

as the averaged value of Xxzx (0) and Xzxx (0) , and hence 
r zxx {0) as the averaged value of r zxx (0) and r xzx (0), for 
all the hexagonal SiC polytypes. Clearly, Table II shows 

(2) 

that Xzxx(0) and hence r zxx (0) are nearly the same for 

(2) 

all the hexagonal SiC polytypes. In contrast, Xzzz(0) and 
hence r 2Z2 (0) increases substantially as one moves from 
2H-SiC to 6H-SiC, suggesting that and r zzz (0) 

are more sensitive to the layer stacking sequence along 
the z-axis. This may be expected since all the electric 
fields are polarized parallel to the z-axis in the latter 
case. 

Our calculated nonvanishing components of x^ (0) f° r 
the four pol ytypes are compared with previous ab initio 
calculations [33, l4ll |42| and available experiments [43l. l44j 
in Table III. Clearly, our results agree rather well with 
the previous ab initio calculations by Chen et al. 14211 . 
Nonetheless, both the results of ours and of chen et aL[43| 
appear to be nearly two-times larger than the previous 
ab initio calculations by Adolph et al. [39] and Rashkeev 
et al. [4l[ . This is surprising since the imaginary parts of 
X^ from both the present and previous [3!| calculations 
are very similar (Figs. 2 and 3). The precise origin of 
these discrepancies is not known at moment. Finally, we 
would consider that the agreement between our results 
and the experimental values is rather good (Table III), 
given the large uncertainties in the available experimen- 
tal data. [H Hi 
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FIG. 3: (Color online) Imaginary parts of x^ 2 \—2lo,lo,uj) of 
SiC in the wurtzite (2H) structure. Previous calculations by 
Adolph et al. [3tJ are also displayed for comparison. 



B. Single graphitic SiC sheet 

In order to understand the nonlinear optical properties 
of the SiC-NTs, we have also calculated the second-order 
nonlinear optical susceptibilities for an isolated honey- 
comb SiC sheet. The isolated SiC sheet is simulated by 
a slab-supercell approach with an intersheet distance of 
about 10 A. The underlying structure was determined 
theoretically by using fc-mesh of 100 x 100 x 1, and the 
theoretical lattice constant is a — 3.069 A. Interestingly, 

(2) (2) (2) 

we find numerically that the x aa b> Xb aa an< ^ Xbbb f° r the 
isolated SiC sheet are nonzero. Here a and b denote the 
two Cartesian coordinates within the SiC layer. Moreo- 

ever > XbL = Xa a \ and xfbb = This is consistent 

with the symmetry consideration, demonstrating again 
that our numerical method and calculations are qualita- 
tively correct. It is because that the isolated SiC sheet 
does not have the spatial inversion symmetry (Dg h ) and 
its symmetry class is D^h (P6m2). Therefore, the iso- 

(2) (2) (2) 

lated SiC sheet has nonzero y n J h , Xbaa anc ^ Xbbb f rom tnc 
symmetry consideration. [381 ] 

We show in Fig. 4 the calculated real and imaginary 
parts as well as the absolute value of the imaginary part 

(2) 

of Xaab(~ w ' w ) °f ^ ne sm gl e SiC sheet. It is clear from 



Fig. 4 that the second harmonic generation (SHG) co- 

efficient X a ab(~ w ' w ) ^ s significant in the entire range 
of the optical photon energy [fkS). Furthermore, for the 

(2) 

photon energy smaller than about 1.2 eV, the X a ab * s 
purely dispersive (i.e., real and lossless) [Fig. 4(a)], sug- 
gesting that the SiC sheet has potential application in 
nonlinear optical devices. Note that since the SiC sheet 
has a theoretical band gap (E g ) of ~2.6 eV, the absorp- 

(2) 

tive (imaginary) part of the Xaab Decomes zer0 below the 
half of the band gap (i.e., ^1.2 eV). The real part of the 

(2) 

Xaab rema ins nearly constant at low photon energies up 
to 0.5 eV, then increases steadily in magnitude as the 
photon energy increases, and finally peaks at the absorp- 
tion edge of ~1.2 eV [Fig. 4(a)] . In the energy range from 

(2) 

1.4 to 3.0 eV, the real part of the X aab becomes positive 
and forms a broad double peak structure. Beyond 3.0 
eV, it becomes negative again and its magnitude grad- 
ually deminishes as the photon energy further increases 
[Fig. 4(a)]. 

(2) 

The absorptive part of the x aa b IS nonzero only above 
~1.2 eV, and looks like a Lorentzian oscillation between 
1.2 and 3.5 eV with one sharp negative peak at ~ 1.4 eV 
and one broad positive peak around 2.8 eV [Fig. 4(a)]. 
It is clear from Eqs. (2) and (3) that the calculated 
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FIG. 4: (Color online) (a) Real and imaginary parts as 
well as (b) the absolute value of the imaginary part of 

Xalb(~ 2w,a>, id) of the isolated SiC sheet. In (c), e"(<^) an( l 
e/Uuj/2) (imaginary part of the dielectric function) from Ref. 
[28 ] are plotted. 



(2) 

X a ab spectra can have pronounced features due to both 
single- and double-frequency resonant terms. To analyze 
the features in the calculated spectra, it is helpful to 
compare the absolute value of x"^ [Fig 4(b)] with the 
absorptive part of the corresponding dielectric function 
e" . Therefore, the calculated e" from our previous pub- 
lication [28| are shown in Fig. 4(c) as a function of both 
uj/2 and id. Clearly, the first sharp peak at ~1.4 eV is 
due to two-photon resonances [cf. e"(u/2)] while in con- 
trast, the second broad peak around 2.8 eV comes from 
the single-photon resonances [cf. e"(w)]- Nevertheless, 
both single- and double-photon resonances involve only 
interband 7r — > 7r* and a — > a* optical transitions for the 
electric field vector E polarized parallel to the SiC layer 
(E || a) M- 

In Table II, the calculated zero frequency linear electro- 
optic coefficient r(0) as well as the corresponding second- 
order nonlinear optical susceptibility x (0, 0, 0) and the 
static refraction index n(0) are listed. The refraction in- 
dex n(0)(= \J e(0)) is derived from the calculated static 
dielectric constant e(0) which has been reported in our 



recent publication 28]. Significantly, the linear electro- 
optical coefficients of the isolated SiC sheet are about 
ten times larger than that of the bulk SiC polytypes (Ta- 
ble II). The static second-order optical susceptibility for 
the SiC sheet is nearly six times larger than the largest 
component of x' 2 ) (0,0,0) of the SiC polytypes. Further- 
more, the zero frequency linear electro-optic coefficient 
r(0) and the corresponding second-order nonlinear opti- 
cal susceptibility x^' (0,0,0) of the SiC sheet are more 
than five times larger than the corresponding values of 
the BN sheet E3. 



C. Second-order optical susceptibility 

We have explicitly calculated the second-order opti- 
cal susceptibility for the zigzag (5,0), (6,0), (8,0), (9,0), 
(12,0), (16,0), (20,0), (24,0), armchair (3,3), (4,4), (5,5), 
(8,8), (12,12), (15,15), andchiral (4,2), (6,2), (8,4), (10,4) 
SiC-NTs. In the case of CNTs, only the chiral nanotubes 
would exhibit second-order nonlinear optical behavior 
with two nonvanishing components of xyz and xzy of 



X 



(2) 



Here z refers to the coordinate along the tube 
axis while x and y denote the two coordinates that are 
perpendicular to the tube axis. As for CNTs, the arm- 
chair SiC-NTs are found not to have any nonzero com- 
ponents of \( 2 \ In contrast, as in the case of the BN- 
NTs (27[, both the zigzag and chiral SiC-NTs are found 
to show second-order nonlinear optical behavior. Specifi- 
cally, all the zigzag SiC-NTs except (5,0), and (9,0), have 
six nonvanishing components of the second-order opti- 
cal susceptibility, namely, xzx, xxz, yyz, zxx, zyy, zzz. 
Nevertheless, these components are not completely in- 

(2) (2) 

dependent of each other. In particular, Xxxz — Xyyz 
Xxzx, and Xzyy 



Xz%- The chiral SiC-NTs have eight 
nonvanishing components of the second-order optical sus- 
ceptibility with two additional nonzero components being 
xyz and yzx. Note that Xyzx — — Xxyz- These numeri- 
cal findings are consistent with the consideration of the 
symmetry of the SiC nanotubes. The point symmetry 
groups of the SiC-NTs [38[ are Ci nv for zigzag (n,0) nan- 
otubes, C-2. n h for armchair (n,ri) nanotubes, and Cjv for 
chiral (n, m) nanotubes where N = 2(n 2 + m 2 + nm)/du 
with dn being the greatest common divisor of 2n + m and 
2m + n. Therefore, these symmetries would dictate [38[ 
that all the components vanish for armchair (n, n) nan- 
otubes, and that nonvanishing components for zigzag 
(n,0) nanotubes are xzx = yzy, xxz = yyz, zxx — zyy, 
ZZZ j cLS well as that nonvanishing components of chiral 
(n, m) nanotubes include all that for zigzag (n,0) nan- 
otubes plus xyz = -yzx. 

We display in Fig. 5 the calculated real and imag- 
inary parts of the second-order optical susceptibility 
•^( 2 )(_ 2id,u>,u>) for the four selected zigzag nanotubes 
[(6,0), (12,0), (20,0), and (24,0)]. As for the single 
SiC sheet, the second harmonic generation coefficients 
(— 2lj, id, id) for the zigzag SiC-NTs are significant in 
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TABLE III: Calculated second-order optical susceptibility X (QL (P m /V) of the four SiC polytypes, together with other 
theoretical [Adolph et al. [39( |. Rashkeev et al. [4ll ]. Chen et al. [43 |] and experimental [Lundquist et al. [43| . Niedermeier 
aZ. results. 





2H 


4H 


6H 


3C 




V (2) Y (2) 
Xzzz-, Xzxx 


(2) (2) 


Y (2) Y (2) 
Xzzz-, Xzxx 


(2) 


This work 


11.5,-17.9 


32.8,-20.4 


38.6,-20.4 


34.2 


Adolph a/. [39] 


4.5, -7.1 


15.5, -9.7 


18.1, -9.8 


17.7 


Rashkeev e£ a/. [41j 


3.6, -6.1 


14.5, -8.9 


17.8, -9.7 


17.5 


Chen et al. [42] 


8.6, -13.2 


23.3, -14.8 


27.6, -15.0 


24.4 


Lundquist et al. [431 






86, -8.6 


17.3 


Nicdcrmcicr et al. [44] 




18±8, -4±2 


24±10, -4±2 






FIG. 5: (Color online) Real (dotted line) and imaginary parts (solid line) of \ ( — 2tt), ui, us) of the zigzag (6,0), (12,0), (20,0) 
and (24,0) SiC nanotubes. 



the entire range of the optical photon energy (huj). In- 
deed, they can be up to ten times larger than that of 
bulk SiC polytypes (see Figs. 1 and 5). Moreover, for the 
photon energy smaller than 2.0 eV, the x' 2 ' is purely dis- 
persive (i.e., real and lossless) (Fig. 5). Interestingly, the 

(2) (2) 

magnitude of Xzzz is the largest (Fig. 5). For Xzzz, the 
electric field of both the incoming and outgoing photons 
is polarized parallel to the tube axis and thus the elec- 
tric depolarization effect would be essentially zero (see 
Ref. |45| and references therein). This may be partic- 
ularly important for nonlinear optical applications. The 
prominant features in each nonzero component of x for 
all the zigzag SiC-NTs except the small-diameter (6,0) 
SiC-NT, look rather similar (Fig. 5). In particular, the 
spectra of the three xzx, zxx and zzz components of 



X (2) for the (20,0) and (24,0) SiC-NTs are nearly identi- 
cal (Fig. 5). We also note that the shape of the spectra of 

(2) (2) 

Xzzz and Xzxx for all the zigzag SiC-NTs look very simi- 
lar, except the difference in sign. For the zigzag SiC-NTs 
with a larger diameter such as (12,0), (20,0) and (24,0), 

the spectrum of Xx'zx is also similar to that of xS'xx (Fig. 
5). The magnitude of the Xxzx and Xzzz components 
decreases slightly as the diameter of the tubes increases 
[e.g., from (6,0) to (24,0)], while the magnitude of the 

(2) 

Xzxx remains more or less independent of the diameter 
(Fig. 5). Furthermore, both the magnitude and shape of 

(2) 

the Xzzz spectrum for the zigzag SiC-NTs with a larger 
diameter, e.g., (20,0) and (24,0), approach to that of the 
single SiC sheet (Fig. 4), as they should. 
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FIG. 6: (Color online) Absolute value of the imaginary part of 
x' 2 ' (— 2uj, lu, (J) (a, c, d) as well as s"(u>) and e"(u/2) (imag- 
inary part of the dielectric function) (b, e) from Ref. [28| of 
the zigzag (12,0) SiC nanotube. 



As for the single SiC sheet, in order to understand the 
features in the calculated x spectra of the zigzag SiC- 
NTs, the absolute values of the imaginary part x ' °f all 
the nonzero components of the (12,0) SiC-NT are plotted 
and compared with the absorptive part of the dielectric 
function e" from our previous publication [2e| in Fig. 6. 
Strikingly, the first prominent peak between 1.0 and 2.0 
cV in the x'zzV spectrum is almost identical to the first 
peak in the e"(oj/2) (see Fig. 6(d) and 6(e)), indicating 
that it is due to two-photon resonances. In contrast, the 
second peak between 2.5 and 3.0 eV in the x"z? spectrum 
is very similar to the first peak in the e"(ui), suggesting 
that it is caused by the single-photon resonances. Never- 
theless, both these single- and double-photon resonances 
involve only interband n —> n* and a — > a* optical tran- 
sitions for the electric field vector E polarized parallel to 
the tube axis (E \\ z) (28|. Fig. 6 also shows that the 
double-peak structure between 1.0 and 2.0 eV in both 
the Xxzx and Xzxx spectra is mainly due to the two- 
photon resonances with E _L z [cf. e"(w/2)] (see Fig. 
6(a)-6(c)), while, in contrast, the second feature in the 



photon energies above 3.0 eV perhaps comes predomi- 
nantly from the single-photon resonances for E _L z [cf. 
e"(w)]. This conclusion is further supported by the fact 
that the magnitude of e"(w) (e"(ui/2)) is only about half 
of that of (e"(w/2)), and concurrently the magni- 

tude of Xxzx and x'zSe is only about half of that of x'zzz 
too. 



In Fig. 7, the calculated real and imaginary parts 
of the second-order optical susceptibility x^ ( — 2w, w, u) 
for all the four chiral nanotubes [(4,2), (6,2), (8,4) and 
(10,4)] are displayed. In general, the spectra of each com- 
ponent of the second-order optical susceptibility for the 
chiral SiC-NTs with a moderate diameter [e.g., (8,4) and 
(10,4)] are similar. For all the chiral SiC-NTs except the 
ultra-small diameter (4,2) SiC-NT, remarkably, the zzz 
component is nearly four times larger than all the other 
nonvanishing components. Another common feature is 
that both the real and imaginary parts of x^ (~2oj, u>, uj) 
show a rather oscillatory behavior, particularly for the 
xyz and xzx components (Fig. 7). The amplitude of 
these oscillatory real and imaginary parts is rather large 
in the photon energies of 1.0~ 4.0 eV for the chiral SiC- 
NTs with a small diameter. It should also be noted that 
the shape and magnitude of the xzx and zxx components 
look very much alike (Fig. 7(b)-7(c)). Finally, both the 

(2) 

magnitude and shape of the Xzzz spectrum for the chiral 
SiC-NTs with a larger diameter, e.g., (10,4), approach to 
that of the single SiC sheet (Fig. 4), as expected. 



Again, in order to understand the structures in the 
calculated spectra of the chiral SiC-NTs, the abso- 
lute values of the imaginary part x"^ °f ai l the nonzero 
components of the (8,4) SiC-NT are plotted and com- 
pared with the absorptive part of the corresponding di- 
electric function e" from our previous publication [28( in 
Fig. 8. Remarkably, the first prominent peak between 
1.0 and 2.0 eV in the x'zzz spectrum is almost identical 
to the first peak in the e"(w/2) (see Fig. 8(d)-8(e)), in- 
dicating that it is due to two-photon resonances. On the 
other hand, the second structure between 2.5 and 3.5 eV 
in the x'zzz spectrum may be correlated with the first 
peak in the £"(w), suggesting that it is caused by the 
single-photon resonances. As for the zigzag SiC-NTs, 
both these single-photon and double-photon resonances 
involve only interband 7r — > 7r* and a — > a* optical tran- 
sitions for the electric field vector E polarized parallel to 
the tube axis (E || z) (28|. Fig. 8 further suggests that 
the feature between 1.0 and 2.0 eV in the spectra of both 
the Xxzx and x'zxx as well as x'xyz may be attributed to 
the two-photon resonances with E _L z [cf. e"(w/2)] (see 
Fig. 8(a)-8(c)), while, in contrast, the second feature in 
the photon energies above 3.0 eV is mainly due to the 
single-photon resonances for E ± z [cf. e"(w)]. 
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Energy (eV) Energy (eV) Energy (eV) Energy (eV) 

FIG. 7: (Color online) Real (dotted line) and imaginary parts (solid line) of x^ 2 '(— 2w, uj, ui) of the chiral (4,2), (6,2), (8,4) and 
(10,4) SiC nanotubes. 



D. Linear electro-optical coefficient 

We list in Table IV the calculated zero frequency linear 
electro-optic coefficient r(0) as well as the corresponding 
second-order nonlinear optical susceptibility x^ 2 K0>0, 0) 
of the SiC-NTs. The r(0) is calculated from the cor- 
responding x^- 2 ' (0)0)0) by using Eq. (7). The refrac- 
tion index n(0)(= ^e(0)) is derived from the calculated 
static dielectric constant e(0) which has been reported in 
our recent publication [28|. Table IV shows that apart 
from the (5,0) and (9,0) SiC-NTs which have no non- 
vanishing r a b c (0), and also the small diameter (6,0) and 
(8,0) SiC-NTs, all the other zigzag SiC-NTs have very 
similar linear electro-optical coefficients, as for the static 
dielectric constant and polarizability [2a ]. Nevertheless, 
r Z zz(0) decreases slightly as the diameter increases, while 
rxzx(0) and r zxx (0) increase slightly as the diameter goes 

(2) 

up. Remarkably, the static Xzzz for the ultra-small diam- 
eter (6,0) is very large (Table IV), being nearly thirteen 
times larger than the largest component of x^ of all the 
four bulk SiC polytypes (Table II). This suggests that 



the SiC-NTs would be useful for applications in the non- 
linear optical devices. Nonetheless, the X zzz decreases 
rapidly as the diamter of the zigzag SiC-NTs increases, 
and, for the large diamter (24,0) SiC-NT, approaches the 

(2) 

value of Xaab °f the single SiC sheet (Table II). The mag- 

(2) 

nitude of the Xaab ^ s about half of the magnitude of the 

X { z% for the (6,0) SiC-NT. We note that X xL and X i% 
are rather similar, and their magnitudes increase as their 
tube diameters increases. Finally, for the larger diameter 

(2) 

zigzag SiC-NTs, the Xzzz is about two times larger than 
Xxzx and X { zL (Table IV). 

The chiral SiC-NTs have two additional nonvanishing 

(2) (2) 

components Xxyz and Xyzx- Nevertheless, the calculated 
static values of Xxyz and Xy 2 zx are negligibly small (i.e., 
within the numerical uncertainty), thereby satisfying the 
requirement by the so-called Kleinman symmetry (40| 

which demands that x2^(0) = Xyzx{Q)- Consequently, 
the corresponding static linear electro-optical coefficients 

(2) (2) 

r xy z(0) and r yzx (0) are zero too. Therefore, Xyzx, Xxyz, 
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TABLE IV: Calculated static refraction index n, second-order optical susceptibility \^ and linear electro-optical coefficient 
r a bc of the zigzag and chiral SiC nanotubes. 



(n z ) 



(2) (2) (2) 
A.XZX-, A.ZXX-, \zzz 

(pm/V) 



' XZX i 1 ZXX j 



(pm/V) 



(5,0) 


2.50 


(3.39) 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 


(6,0) 


2.53 


(3.31) 


77.3, 45.5, -495.5 


-2.19, -1.29, 8.21 


(8,0) 


2.56 


(3.23) 


87.6, 68.0, -323.1 


-2.56, -1.98, 5.92 


(9,0) 


2.58 


(3.23) 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 


(12,0) 


2.60 


(3.21) 


98.9, 88.1, -255.4 


-2.84, -2.53, 4.82 


(16,0) 


2.62 


(3.21) 


104.7, 96.9, -241.9 


-2.96, -2.74, 4.57 


(20,0) 


2.68 


(3.27) 


125.0, 119.6, -279.2 


-3.25, -3.11, 4.85 


(24,0) 


2.63 


(3.19) 


107.1, 102.1, -227.7 


-3.05, -2.91, 4.39 


(4,2) 


2.74 


(3.37) 


23.9, 4.46, -410.4 


-0.56, -0.11, 6.36 


(6,2) 


2.56 


(3.16) 


58.2, 37.0, -236.6 


-1.78, -1.13, 4.75 


(8,4) 


2.62 


(3.17) 


50.3, 43.6, -138.6 


-1.46, -1.26, 2.75 


(10,4) 


2.62 


(3.17) 


65.3, 59.3, -164.3 


-1.89, -1.72, 3.25 



r xyz(0) and r yzx (0) for the chiral SiC-NTs are not listed 
in Table IV. As for the nonvanishing static components, 
r xzx (0) and r zxx (0) for the (6,2), (8,4) and (10,4) SiC- 
NTs are quite close. xizz(O) and r zzz (0) decrease as the 

(2) 

tube diameter increases, while, in contrast, Xzxx(0) and 
r zxx (0) increase with the tube diameter. r xzx (0) and 
r zxx (0) for the (4,2) is somewhat smaller than that of 
the other chiral SiC-NTs. 

The zero frequency second-order nonlinear optical sus- 
ceptibility x< 2 )(0, 0, 0) of the SiC-NTs (Table IV) are gen- 
erally a few times larger than the counterparts of the cor- 
responding BN-NTs [13 . In particular, the xfzz (0,0,0) 
of the (6,0) SiC-NT is about thirteen times larger than 
that of the (6,0) BN-NT [27j]. On the other hand, the 
xi 2 2(0, 0, 0) of the (24,0) SiC-NT is only about six times 
larger than that of the (24,0) BN-NT [27j]. This suggests 
that compared with BN-NTs, SiC-NTs would be a better 
nonlinear optical material. However, the low frequency 
linear electro-optic coefficient r(0) of the SiC-NTs (Table 
IV) are only slightly larger than the counterparts of the 
corresponding BN-NTs [27j. For example, the r zzz (0) of 
the (6,0) and (24,0) SiC-NTs is, respectively, 1.7 and 1.3 
times larger than that of the (6,0) and (24,0) BN-NTs. 
This is because the refraction index n of the SiC-NTs 
(Table IV) are also larger than that of the BN-NTs H3 
[seeEq. (7)]. 



(6,0), (8,0), (9,0), (12,0), (16,0), (20,0), (24,0)], armchair 
[(4,4), (5,5), (8,8), (12,12), (15,15)], and chiral [(4,2), 
(6,2), (8,4), (10,4)] nanotubes have been calculated. For 
comparison, the second-order nonlinear optical proper- 
ties of bulk SiC polytypes (2H-, 4H-, 6H- and 3C-SiC) 
and the single graphitic SiC sheet have also been cal- 

(2) (2) 

culated. Interestingly, we find that the X a ab> Xbaa ano - 

(2) 

Xbbb f° r * ne isolated SiC sheet are large and generally 
several times larger than that of the SiC polytypes. Un- 
like carbon nanotubes, both the chiral and zigzag SiC- 
NTs have pronounced second-harmonic generation and 
linear electro-optical coefficients which are comparable 
to or even larger than that of the single SiC sheet. The 
prominant structures in the spectra of x^ ( — 2u>j u>, uj) of 
the SiC-NTs have been related to the features in the cor- 
responding linear optical dielectric function e{u>) in terms 
of single-photon and double-photon resonances. We also 

(2) 

find that the X a b c ano - r<lbc coefficients of the SiC-NTs are 
up to thirtcen-times larger than the counterparts of the 
corresponding BN-NTs. Therefore, the SiC-NTs would 
be promising nonlinear optical materials for applications 
in, e.g., second-harmonic generation, sum frequency gen- 
eration and electro-optical switches. We hope that this 
work would stimulate experimental investigations into 
the second-order nonlinear optical properties of the SiC- 
NTs. 



IV. SUMMARY 
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FIG. 8: (Color online) Absolute value of the imaginary part 
°f X*' 2 ' 1 (~2LJ,tJ,tj) (a, c, d) and e"(u) and e"(uj/2) (unaginary 
part of the dielectric function) (b, e) from Ref. [28|] of the 
zigzag (8,4) SiC nanotube. 
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